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Abstract 

Lattice structures due to their low weight and high 
performance as structural elements have been widely used 
in different aerospace applications. In this paper, effective 
parameters on the design of anisogrid lattice conical shells 
are investigated. First, filament winding patterns regarding 
the desired axial strength is decided. Then, regarding 
geometrical relations, effective parameters in order to form 
the anisogrid cell are identified. Distance of circular ribs 
from each other has an important role in determining the 
conical lattice structure and eventually in deriving the 
stiffness matrix. Finally, considering the relations, finite 
element analysis model of the lattice conical structure has 
been performed by ABAQUS software and buckling analysis 
under axial loading is done. In deriving the strength results, 
the verifications of references of this paper and the classic 
theory have been used. 
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Introduction 

Composite lattice structures are one of the most 
complex and newest structures, which are designed 
and produced in both isogrid and anisogrid forms. In 
the past decade, researches done on polymeric matrix 
fiber composite lattice structures have become one of 
the scientific and widely used centers of attention, and 
these structures have a widespread application in 
different aerospace structures. These structures as 
lattice lamina are used lonely or with internal or 
external shells. Lattice lamina in these structures 
includes different systems of ribs, which are produced 
from continuous fibers, using automatic filament 
winding method. In these structures under arbitrary 
loading fibers of the lattice are under tension or 
compression, which is a quite desirable state for 
composites, and is one of the major advantages of 
these structures compared to the conventional 



structures for composites. Major characteristics of 
these composites in the face of axial compression and 
bending moment loadings, are different failures 
occurring due to total buckling of the structure, failure 
due to maximum stress in helical ribs, and failure due 
to local buckling of the helical ribs in the lattice. 
Buckling load related to each of the failure states 
depends upon geometric parameters and total 
dimensions of the structure. In the analysis of these 
structures it is supposed that elements of the lattice are 
two force members, behaving as an orthotropic body, 
and the relations of orthotropic structures are used 
here for the analysis of the total buckling of these 
structures. Major application of these structures is in 
aerospace industry, which widely uses different forms 
of them, including: curved shaped and annular shells 
(especially conical shells). For the purpose of 
improving the mechanical properties of the structure, 
and its weight optimization, composite materials have 
been used. In the midst of 1960's the finite element 
method was developed for the numerical analysis of 
these structures. Generally, these structures are 
manufactured by advanced forming of fibers. Also, the 
ribs use in these structures can move in 2 up to 4 
directions. Design, analysis and manufacture of these 
structures are investigated in reference [1]. More 
details of their analysis are included in references [2, 3]. 
In references [4, 5] buckling load analysis of grid 
stiffened cylindrical shells is performed. Also, the 
results of the experimental test are compared and 
verified in comparison with the analytical results. Kim 
[6, 7] has studied the fabrication and testing of thin 
composite isogrid stiffened panels and composite 
isogrid stiffened cylinders, so as to investigate their 
buckling behavior. In reference [8] design and 
optimization of laminated conical shells for buckling 
and maximum buckling loads has been performed. In 
this study, optimization has been performed in two 
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states of causing the maximum buckling load at 
certain weight, and causing the minimum weight 
under a constant critical load. In reference [9] critical 
buckling load is derived from the solution of the 
governing nonlinear partial differential equation with 
different coefficients. Reference [10] has investigated 
the optimized stacking sequence design, so as to 
achieve the maximum fundamental eigenfrequency in 
conical shells. In reference [11] buckling of the 
cylindrical lattice structure has been done, using finite 
element method. Also, in reference [12] finite element 
analysis for buckling of the conical lattice shell is 
studied. In this paper, a numerical code is developed 
for the conical lattice structure, which determines the 
basic parameters, buckling under critical axial loading, 
torsion moment, and bending moment. In this paper, 
design of the conical lattice structure with the 
anisogrid cell is investigated. First, the governing 
differential equations of the anisogrid cell of the 
conical lattice structure are derived, regarding the 
filament winding technique used. Then variations of 
the principal parameters of the design are investigated 
in relation to the increase of the number of helical ribs. 
Later, by creating a finite element model of the conical 
lattice structure with the anisogrid cell, buckling 
analysis of the structure under axial loading is 
performed. Also, the variations of the critical buckling 
load, relative to increase of thickness, width, and the 
distance between ribs is considered. 

Governing equations of design 

Lattice structures are widely used in different 
industries. These structures are highly strong 
against the destruction caused by impact, 
delamination, and craze propagation in the 
structure. One of their advantages over other 
structures is reduction of weight and manufacture 
time. Lattice structure is made of a number of 
helical and circular ribs. These structures can have 
up different forms, including isogrid and anisogrid. 
Generally, the main purpose of using lattice 
structures is optimal usage of longitudinal 
properties of the composite materials used. Also, a 
schematic of the conical lattice shell used as a 
satellite carrier adaptor, at the lower section of the 
satellite is shown in figure 1. 

All the parameters are determined on the basis of the 
location of helical and circular ribs. Some of these 
parameters depend on the angel of the helical ribs. 
Also, some of these parameters as the height of the 



structure increases have variations with a certain rate. 
Because of the complexity of the geometry of conical 
lattice shells, in order to design, one first must divide 
the effective design parameters of the structure in two 
dependent and independent categories, and then 
perform the complete investigation procedure of 
design. 



Composite Shroud 

X 




FIGURE 1 a SCHEMATIC OF THE SATELLITE CARRIER 
ADAPTOR [13] 

Independent parameters 

The number of circular and helical ribs (n c , n h) are 
independent parameters in the design of conical lattice 
structure. 

Dependent parameters 

Dependent parameters include: Ai// < (p < X < a c 
< «/i • 

A schematic of the parameters stated here in the 
design of annular lattice structure, is shown in figure 2. 




FIGURE 2 GEOMETRIC PARAMETERS OF THE LATTICE 
STRUCTURE 

Since the design of these structures is intended to 
stand axial compression loadings, therefore the most 
important problem regarding these structures is the 
analysis of buckling and elastic stability. Considering 
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these facts, design of the lattice structure and 
eventually manufacture of it, must be based on a state 
that shows the maximum strength against buckling. 
Amongst different filament winding patterns, geodesic 
pattern is the best pattern of filament winding for 
lattice structures, carrying axial compression loadings. 
In order to derive the governing differential equations 
of the structure, first by developing a section of the 
cone, the geometric equations, which have a direct 
effect on the stiffness matrix, are achieved (figure 3). 



assumption that the circular ribs are located 
between the intersection of the helical ribs. 
Considering the geometry of figure 3, the following 
equation is achieved: 



y = -?-{n -1) 



(4) 



Geodesic angle at the smaller cross-section ( (p Q ) is 
computed as follows: 




<p = tan 



p F sin y 
p F co&Y-p 



(5) 




Pn 




V 





FIGURE 3 DEVELOPED SECTION OF THE CONICAL 
LATTICE SHELL 

Using geodesic relations, the following equation 
for each point of the lattice structure exists [6]: 



p sin cp = p o sin q> o = C o = constant 



(1) 



Wherein 



P 



and r is the radius of the 



sin a 

cross section of the cone which is different at each 
point, and OL is the angle between the line passing 
through the apex of the cone and the axis of 
symmetry of the cone. Taking derivation of 
equation (1) gives: 



dp 
P 



dq> 
tan cp 



(2) 



Where Ay/ is the angle between the two helical 
ribs of figure 3, and is found here: 



Ay/ = 



2tt sin a 



(3) 



Geometrical design and finite element analysis of 
the conical lattice structure is performed with the 



And also, the geodesic angle at the bigger cross-section 
( ^> ) is equal to: 



<Pf = <P* ~ 1 



(6) 



Variations of p versus (p are achieved by this 
equation: 



dp = -C„ C ° dcp 
sin (p 



(7) 



Vertical space between circular and helical ribs for a 
cell at each row is derived from the following 
equations: 



a h - 2a c sin cp 



(8) 



f A+1 A y/dp 
hh 

r +1 2d(p 



(9) 



In equation (9), i indicates the number of ribs. The 
minus sign in equation (9) shows that the axis of the 
system of coordinate system is located on the smaller 
cross-section. Stiffness properties of one cell, which 
has a repeating pattern, is the representative of the 
total repeating section. The orthotropic properties 
achieved for annular bodies, is along the axial 
direction. For example, for the cylindrical structure, 

— A 

stiffness along the axial direction is equal to E x — , 

L 

where A = 2tiRH is the cross-section of the structure. 
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(10) 



5 = smcp ,c — cos (p 

E h and £ are the moduli of elasticity of helical and 

circular ribs, respectively. Equation (10) is directly 
derived from theory, formulation and assumptions, 
which are related to the stiffness of the layer and the 
properties of fibers. Also, properties of the equal 
stiffness along the axial direction are derived from 
equation (11), [6]: 



E = 



H 



q n q 22 



2 \ 

In 



In 



(11) 



^■used in equation (11) are the components of the 
stiffness matrix Q , and as it is seen, they depend upon 
the angle of helical ribs, width of the ribs, and the 
distance between helical and circular ribs. The critical 
buckling axial load for conical shell is given by [14]: 



cylco 



2nE X 2 
-cos a 



V3d-v 2 ) 



(12) 



According to equation (12), the critical buckling axial 
load for conical shells is the same as that of cylindrical 
shells, with this difference that the angle for half-apex 
of the cone is also effective in the critical load. By 
comparison of the experimental, analytical and finite 
element analysis, it is found out that in order for the 
analytical result to have an acceptable answer, a 
buckling correction factor of C must be multiplied 
with it. Eventually, by applying the buckling 
correction factor, the critical buckling load of the 
conical shell is equal to: 



P = C 



2nEf 2 

^cos a (13) 



a/3(1-v 2 ) 



According to reference [14], the magnitude of 
parameter C depends on the angle of the cone's 
inclination. The correction factor of C for a cone, 
having the apex angle of 10 up to 75 degrees, is equal 
to 0.33. At a range outside these angles, C is derived 
according to experimental test. 



Finite element model 

Regarding the fact that these structures cannot be 
transformed into smaller samples, in order to 
determine the mechanical specification of the structure, 
one can only use the possible solution of experimental 
test. Numerical analysis of the lattice structure is 
performed, using Timoshenko's beam element, and the 
shell's element, as it is seen in figure 4. Generally, 
considering the situation that major loads applied to 
the lattice structure are axial loads, therefore buckling 
is the most important parameter in strength evaluation 
of these structures. Choosing of the beam element is 
described in reference [11, 12]. Circular and 
circumferential ribs of the lattice structure are 
completely bonded to each other at the intersection 
points. The distance between circular and helical ribs 
from each other is achieved, considering the governing 
relations. 




FIGURE 4 A SCHEMATIC OF THE BEAM ELEMENT USED IN 
FINITE ELEMENT ANALYSIS 

In order to design the winding pattern of the fibers, 
it is necessary to derive the geodesic angles of the 
helical ribs in big and small cross-sections, and the 
height location of the conical structure should be 
known. For the design of conical lattice structure, 
and regarding equations (6), (7), and (10) the total 
dimensions including height, bigger and smaller 
cross-sections, and the apex angle must be 
determined. By this explanation, big and small 
diameters of the cone are 2500 mm and 1250 mm, 
respectively, and the apex angle is 34.7°. Now, 
regarding the geometric specification stated, 
geodesic angle variations at different height 
locations sand different number of circular ribs is 
derived. Considering the model at the previous 
section, geodesic angle variations at the bigger and 
smaller cross-sections with the increase of the 
height of the cone, and variations of the number of 
circular ribs, are shown in figure 5 and 6, 
respectively. As it is seen, with the increase of the 
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height of the cone, geodesic angles are non-linearly 
decreasing. Also, with the increase of the number 
of circumferential ribs at the same geometric 
conditions, geodesic angles at each section are 
reduced. Results show that by increasing the 
number of circumferential ribs, variations of the 
geodesic angle at the bigger section are reducing 
with a more inclination, compared to the smaller 
section. Regarding the total geometric 

specifications of the conical structure, geodesic 
angle can be achieved at the bigger and smaller 
cross-sections for a certain number of helical and 
circular ribs from the diagrams of figure 5 and 6. 
For example, for the conical lattice structure with 

10 circular ribs, and the ratio of — geodesic 

P. 

angle at the bigger and smaller cross-sections are 
15.98 and 33.44, respectively. Also, variations of X 
with the increase of the number of helical ribs, is 
shown in figure 7. The finite element model with 
total specifications and geodesic angle is shown in 
figure 8. 




a J ■ 1 1 % 



FIGURE 5 VARIATION OF <p o WITH THE INCREASE OF B±- 
(a = 34.7°) 




FIGURE 6.VARIATION OF (p F WITH INCREASE OF tllL 

Po 

(a = 34.7°) 





2J . 




22 ■ 


-*- 


20 ■ 




18 ■ 










16 ■ 


X 






1-1 ■ 




12 




in ■ 



40 5(1 60 "0 80 

FIGURE 7 VARIATION OF I WITH THE INCREASE OF THE 
NUMBER OF HELICAL RIBS 




FIGURE 8 FINITE ELEMENT MODEL OF THE LATTICE 
CONICAL SHELL, MADE BY ABAQUS SOFTWARE 

Verification 

Finite element model of the lattice conical shell which 
has undergone buckling analysis under axial and 
bending loadings in reference [12] is used here and 
buckling analysis is performed on it, using ABAQUS 
software and Timoshenko's beam element. Geometric 
specifications and the materials used in the finite 
element analysis are chosen according to reference [12]. 
A schematic of the finite element analysis of the lattice 
conical shell with the defined geometric specifications 
is shown in figure 9. Critical loads here, are in Newton, 
and for better display of the deformation of the 
structure, the aspect ratio of 40 is applied. 

Results presented in table 1 show that the results 
achieved from the numerical finite element 
solution and the analytical results are of high 
accuracy. Regarding the verification performed in 
the following on the lattice conical and cylindrical 
structures with certain geometrical specifications, 
and the equations stated above, weight 
optimization of the structure is performed. 
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TABLE 3 PROPERTIES OF THE MATERIALS USED IN THE 
LATTICE CONICAL SHELL 



(GPa) 


(MPa) 


(Vs) 

/ m 


E c 

(GPa) 


m c 

(V 3 ) 

/ m 


80 


350 


1450 


64 


1410 



According to the diagram of figure 10, analytical 

a c 

results show that by increasing the — ratio, critical 

a h 

buckling load increases. 



■Step: Stcp-1 

Mqd- 1 ■ Eig*r>V*m » 15.039 

PnimurvVcir; U, Muwiiludu 



FIGURE 9 DERIVATION OF THE CRITICAL BUCKLING 
LOADING, USING FINITE ELEMENT METHOD (H=B=4MM, 

a = 5°) 



TABLE1 COMPARISON OF THE RESULTS OF REFERENCE 
[12] WITH THE RESULTS OF ANALYTICAL ANALYSIS, AND 
THE NUMERICAL ANALYSIS PRESENTED HERE 



h=8mm 
b=2mm 


h=4mm 
b=4mm 


h=2mm 
b=8mm 




44.334 


34.416 


13.863 


Reference 
[12] 


49.112 


33.831 


15.03 


ABAQUS 


47.997 


36.216 


15.06 


Analytical 



Results and discussion 

Critical buckling load in structures depends upon 
stacking sequence of the ribs, angle of helical ribs, 
cross-section of the dribs, and etc. Later in this 
paper, results of variations of axial loading under 
clamped support conditions are investigated. 
Geometric specifications of the lattice conical shell 
are summarized in table 2 with regards to figures 5 
and 6. Also, the properties of materials used in the 
lattice conical shell are included in table 2. 

TABLE 2 GEOMETRIC PARAMETERS OF THE LATTICE 
CONICAL SHELL STUDIED 





n h 


a f 




H 

(mm) 


K 

(mm) 


h 

(mm) 


10 


53 


15.98 


33.44 


18 


4 


5.75 
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FIGURE 10 VARIATION OF BUCKLING LOAD, VERSUS 

a,. 



INCREASE OF- 



a. 



Also, with the increase of — ratio, critical buckling 

K 

load decreases (figure 11). Also, variations of the 
critical buckling load, by increasing the — ratio, non- 

K 

linearly decreases (figure 12). In figure 13, variation of 
the critical buckling load versus increase of thickness, 
has been shown for two states of analytical and finite 
element methods. Results indicate that, by increasing 
the thickness of the ribs of lattice shell, the critical 
buckling load of the structure, non-linearly increases. 



LODE -09 -I 
8.00E+Q8 - 
g 6.0DE+0S ■ 

Ffl 

ft 4.00E+08 ■ 
2.O0E+O8 ■ 
O.OOE+00 • 




a c /b c 



FIGURE 11 VARIATION OF CRITICAL BUCKLING LOAD, 

a c 

VERSUS INCREASE OF RATIO 

b. 
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FIGURE 12 VARIATION OF CRITICAL BUCKLING LOAD, 



VERSUS INCREASE OF 



RATIO 



'FEW 



Analytical 




QQDE+00 



5 10 15 :o 

Thickness {mm) 

FIGURE 13 VARIATION OF CRITICAL BUCKLING LOAD, 
VERSUS INCREASE OF THICKNESS OF CIRCULAR AND 
HELICAL RIBS 

Regarding the diagram of figure 13, it is realized that 
with the increase of the thickness of the rib 
(approximately 20mm), the beam element does not 
possess high accuracy in the convergence of the results 
of analytical and numerical analysis. 




10mm 



*" rr-ryffir"'~~. 



5mm 





20mm 



15mm 



FIGURE 14.A SCHEMATIC OF THE FIRST MODE OF BUCKLING 
FOR THE BUCKLING OF LATTICE CONICAL SHELL UNDER 
AXIAL LOADING, VERSUS THICKNESS VARIATION OF THE 
RIB (WIDTH=5.75MM) 



Also, with the increase of the width of helical and 
circular ribs, critical buckling load increases non- 
linearly (figure 15). 



■ FEM 



-Analytical 




5 10 15 

Width (mm) 

FIGURE 15 VARIATION OF CRITICAL BUCKLING LOAD, 
VERSUS INCREASE OF WIDTH OF CIRCULAR AND 
HELICAL RIBS 

Regarding the results presented, one can derive the 
critical axial load with a good accuracy, using 
numerical finite element method and analytical 
method. A schematic of the first up to sixth modes of 
buckling of the lattice conical structure, under axial 
loading is presented in figure 16. 





20mm 

FIGURE 16 FIRST MODE OF BUCKLING FOR THE 
LATTICE CONICAL SHELL UNDER AXIAL LOADING, 
VERSUS WIDTH VARIATION OF THE RIB 
(THICKNESS=18MM) 



Final conclusion 

Design of a lattice conical shell, requires considering 
different parameters. These parameters are 
independently chosen by the designer, or are 
dependent. For buckling analysis of the lattice conical 
structure under axial loading, the equal stiffness of the 
shell in the axial, circumferential, and through the 
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thickness directions should be determined. Stiffness 
matrix of the lattice structure, besides the mechanical 
properties of helical and circular directions, depends 
upon the geometry of the cells of the lattice structure, 
and also the angle of the helical ribs at different points 
of the height. In buckling analysis of the lattice conical 
shell under axial loading, the following results are 
achieved: 

a c 

■ By increasing the — ratio, the critical buckling 

a h 

load increases. 

a c 

■ By increasing the — ratio, the critical buckling 

K 

load decreases. 

a h 

■ By increasing the — ratio, the critical buckling 

K 

load decreases. 

■ By increasing the thickness of the circular and 
helical ribs, the critical buckling load increases 
non-linearly. 

■ By increasing the width of the circular and helical 
ribs, the critical axial loading increases. 
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